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I. INTRODUCTION 



Chiral effective field theory (EFT) provides a systematic and model-independent framework to analyze low-energy 
hadronic processes in harmony with the spontaneously broken approximate chiral symmetry of QCD. This approach 
has been successfully applied to the derivation of the nuclear forces and, more recently, also hyperon-nucleon and 
hyperon-hyperon interactions, see Refs. [l|,|2|, |3[ for review articles. Exchange vector and axial currents in nuclei have 
also been studied in the framework of chiral EFT. In their pioneering work, Park, Min and Rho applied heavy-baryon 
, chiral perturbation theory to derive exchange axial Q and vector Q currents at the one-loop level for small values 
of the photon momentum focusing, in particular, on the axial-charge and magnetic moment operators, respectively. 
These calculations were carried out employing time-ordered perturbation theory to extract non-iterative contributions 



to the amplitude. The resulting exchange vector currents were applied within a hybrid approach to analyze magnetic 
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moments and radiative capture cross sections of thermal neutrons on light nuclei [a, B 0> HI as well as some polarization 
observables in radiative neutron capture on the proton Q. For applications to various electroweak few-nucleon 
reactions of astrophysical interest see [l(| EH, EH EH, E3 • 

Deuteron electro mag netic properties [H, EH E3, EM EH, Compton scattering on the deuteron [2(], HH and, more 
' recentl y, o n 3 He [22j | as well as pion electro- and photoproduction and the corresponding capture reactions [23l . l24l . 
HI, l26Ll27l |28j have also been addressed in the framework of chiral EFT. However, to the best of our knowledge, no 
applications have so far been performed to electron and photon inelastic few-nucleon reactions with the momentum 
transfer of the order M r where a lot of experimental data are available, see [2!| for a recent review article on the 
theoretical achievements in this field based on conventional framework. Recent progress in the accurate description 
of the two- [3(| HH and more-nucleon systems [13] within the chiral EFT, see also Q and references therein, gives 
a strong motivation to apply this framework to the abovementioned processes. This requires the knowledge of the 
consistent electromagnetic exchange current operator for non-vanishing values of the photon momentum. While the 
leading two-nucleon contributions to the exchange current arise from one-pion exchange and are well known, the 
corrections at the one- loop level have not yet been completely worked out. An important step in this direction was 
done recently by Pastore et al. [H, HH who considered the electromagnetic two-body current density at the leading 
one-loop order based on time-ordered perturbation theory. In the present work, we calculate the leading two-pion 
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exchange two-nucleon four-current operator in chiral EFT based on the method of unitary transformation which 
we used to derive nuclear forces in Refs. [H, [H, H3, [H, [39j|. Our work provides an important check of the results 
presented in Refs. [H, H3] but also differs from these works in several important respects. First, as already pointed 
out, we use a completely different method to compute the current operator. Secondly, we also give results for the 
exchange charge density which, to the best of our knowledge, have not yet been published before. Finally, we evaluate 
analytically all loop integrals to obtain a representation in momentum space in terms of the standard loop functions 
and the three-point functions. The latter are reduced to a form which can be easily treated numerically. We also 
succeeded to analytically carry out the Fourier transformation for all contributions leading to an extremely compact 
representation of the current and charge densities in coordinate space. Notice that contrary to [H, [Hj], we do not 
treat the A(1232) isobars as explicit degrees of freedom in this work. 

Our manuscript is organized as follows. In section [IT] we provide a short summary of the method of unitary trans- 
formation, explain the adopted power counting scheme, list all relevant terms in the effective chiral Lagrangian and 
present our results for the exchange current and charge densities in momentum space. The expressions in coordinate 
space are given in section IIII1 The results of our work are summarized in section IIVI The formal operator struc- 
ture of the effective electromagnetic current can be found in appendix \X\ while appendix [B] contains the complete 
momentum-space representation of the current. Evaluation of the three-point functions entering this representation 
is detailed in appendix [Cl In appendix [D] we verify that the obtained expressions for the current fulfill the continuity 
equation. Finally, appendix [E] collects the formulae needed to obtain the coordinate-space representation given in 
section IIIII 



II. NUCLEAR CURRENTS USING THE METHOD OF UNITARY TRANSFORMATION 



We begin with a brief reminder about the method of unitary transformation. Consider the time-independent 
Schrodinger equation for interacting pions and nucleons in the absence of electromagnetic sources 

(ffo + H/)!*} =E\9), (2.1) 

where |^) denotes an eigenstate of the Hamiltonian H with the eigenvalue E. Let T) (A) be projection operators onto 
the purely nucleonic (the remaining) part of the Fock space satisfying rj 2 = r), X 2 = A, r/X = Xrj = and A + r) = 1. 
To describe the dynamics of few- and many-nucleon systems below the pion production threshold it is advantageous 
to project Eq. (|2.1|1 onto the 77-subspace of the full Fock space. The resulting effective equation can then be solved 
using the standard methods of few- or many-body physics. The decoupling of the rj- and A-subspaces can be achieved 
via an appropriately chosen unitary transformation [40L l4lj 

jg, (2.2, 

see Ref. [42| for an alternative approach. Following Okubo [4l|, the unitary operator U can be parametrized as 

_(r,(l + AU)-V= -At(l + AAt)-V2 \ 
U ~ { A(l + AU)-V2 A(l + yL4t)-i/2 ) > V- 6 ) 

in terms of the operator A = XArj which has to satisfy the decoupling equation 

A (H - [A, H] - AHA) rj = (2.4) 

in order for the transformed Hamiltonian H to be of block-diagonal form. The effective ry-space potential V can be 
expressed in terms of the operator A via: 



V = V (H- H 



(1 + A^A)-^ 2 {H + A^H + HA + A*HA)(1 + A t A)- 1 ' 2 - H 



(2.5) 



The unitary transformation U and the effective potential V can be calculated perturbatively based on the most general 
effective chiral Lagrangian utilizing the chiral power counting. In Ref. [39], a convenient formulation of the power 
counting has been presented. The low-momentum dimension v of the effective potential, V,V~ 0(Q/A) u with Q and 
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A refering to the soft and hard scales of the order of the pion and p-meson masses, respectively, is given (modulo the 
normalization constant —2) by the overall inverse mass dimension of the coupling constants entering the expression 
for V: 



^n t + Pl ~A. (2.6) 



Here, where Vi is the number of vertices of type i while di, rij andft refer to the number of derivatives or M^-insertions, 
nucleon field and pion field operators, respectively. Further, «j is simply the canonical field dimension of a vertex of 
type i (up to the additional constant —4). Writing the effective chiral Hamiltonian H as 



H = J2 H(K) > ( 2 - 7 ) 



the operator A can be calculated recursively, 



a — 1 ' 



-1 a — 1 a — 2 a — j — 1 



-|- y^ jjW^4( Q; -') _ y~^ ^(»-0g(*) _ y~^ y^ ^WjyO')^"-*-.?) 

i—1 i—1 i—1 j — 1 



77. (2.8) 



Here, (-E\) refers to the free energy of nucleons (nucleons and pions) in the state rj (A). The expressions for 
the unitary operator and the effective potential then follow immediately by substituting Eqs. (|2.7p and (|2.8p into 
Eq. (|2.5[) . It is important to emphasize that Eg. (|2 . 3[) does not provide the most general parametrization of the 
unitary operator. Moreover, as found in Ref. [39(, the subleading contributions to the three-nucleon force obtained 
using the parametrization in Eq. (|2.3p cannot be renormalized. To restore renormalizability at the level of the 
Hamilton operator additional unitary transformation U' in the 77-subspace of the Fock space had to be employed, 
rjU'rjU rj = rjU rjU'r] — 77, whose explicit form at lowest non-trivial order is given in that work. 

It is, in principle, straightforward to extend this formalism to low-energy electromagnetic reactions such as e.g. electron 
scattering off light nuclei, see [43|, HU, 45] for some early applications of the method of unitary transformation to the 
derivation of the exchange currents. Here and in what follows, we restrict ourselves to the one-photon-exchange 
approximation to the scattering amplitude. The effective nuclear current operator 77 (2)77 acting in the 77-space is 
then defined according to 

{*,\Ji>(x)\*i) = (folvU'^UU^UvU'r)^) = (tPflnJ^xWfr) , (2.9) 

where ?7|<^j,/) = njU rjU^l^ij) denote the transformed states and we have omitted the components X\<j>if) which 
is justified as long as one stays below the pion production threshold. In the above expression, J ll (x) denotes the 
hadronic current density which enters the effective Lagrangian C^Ny describing the interaction of pions and nucleons 
with an external electromagnetic field A^ and is given by 

"W-^fe-T^ 1 - <2 - 10) 

Notice that contrary to the Hamilton operator, the unitarily transformed current does, in general, not have the 
block-diagonal form, i.e. r/W J^(x)UX 7^ 0. Again, it is important to realize that the above definition of rjJ^ s (x)ri 
does not fully incorporate the freedom in the choice of unitary transformations. Thus, one might expect that this 
formulation yields the effective current operator which is not renormalizable by a redefinition of the low-energy 
constants (LECs) entering the underlying Lagrangian. Indeed, renormalizability of the effective current operator 
implies highly non-trivial constraints in the case of one-pion exchange contributions at the one-loop level since all 
/3-functions of the LECs U from [49|, [5(| and di from C^n f5~Ll |52| are fixed. We have verified that the ultraviolet 
divergences entering the expressions for the one-pion exchange contributions using the formulation based on the A 11 - 
independent unitary transformation as described above can indeed not be completely removed by the redefinition of 
the corresponding LECs. Thus, a more general parametrization of the unitary transformation is required in order 
to restore renormalizability of the nuclear current. This can be achieved if one allows for the unitary operator to 
depend explicitly on the external electromagnetic field, U(A^ L ). The operator U{A^) then has to be chosen in such 
a way that the transformed Hamiltonian W (Af 1 ) H^n-v U (A^ ) is block-diagonal (with respect to the rj- and A-spaces) 
and coincides with the one given in Ref. [39| when the external electromagnetic field is switched off. The effective 
nuclear current operator rjj^ s {x)rj in this more general formulation receives additional contributions which are not 
included in Eq. (|2.9p and result from ^-dependent pieces of U(A^) in the expression (A fJ ')H 7! -N 1 U(A ,J ') whose 
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form is determined by renormalizability of the resulting nuclear current operator. These additional terms in r]J^ s (x)rj 
are found to have no effect on the two-pion exchange current and will be discussed in detail in a separate publication 
devoted to the one-pion exchange contributions. Finally, we emphasize that the power counting employed in the 
present work implies the following restrictions on the photon momentum k in the two-nucleon rest frame 

|fc| ~ O (M„) , k ~O ( J « Af„ , (2.11) 

where and m refer to the pion and nucleon masses, respectively. For the kinematics with fc° ~ O (J14), one will 
have to systematically keep track of the new soft momentum scale y/M^m. This goes beyond the scope of the present 
work. 

For the calculation of the leading two-pion exchange two-nuclear current operator in the present work we only need 
the leading pion and pion-nucleon terms in the effective Lagrangian 

= !ftt[DpUiy*tf + M$(U + tf)] , 

4w = N\ivD + g A S-u)N, (2.12) 

where the superscript i in £w denotes the number of derivatives and/or pion mass insertions. Here, F v (qa) is the 
pion decay constant (the nucleon axial- vector coupling) , N represents a nucleon field in the heavy-baryon formulation 
and 5 M = Yl^a^vV" is the Pauli-Lubanski spin vector which reduces to = (0, for — (1, 0, 0, 0). At the order 
we are working and for the contributions to the current operator considered in the present work, all LECs entering 
Eq. (|2.12|) should be taken at their physical values. Further, the SU(2) matrix U — u 2 collects the pion fields and 
various covariant derivatives are defined according to 

D„U = duU-irpU + iUlp, 



= i [u^(dfj, - ir^)u ~ u(<9 M - i/ M )u t ] 



fl M + r^-i«W]jV with T^ = ^[u\d^-ir^)u + u{d^-il^} . (2.13) 



To describe the coupling to an external electromagnetic field, the left- and right-handed currents and and the 



isoscalar current v^f 1 have to be chosen as 



I* = , «W = e A , , (2.14) 



where e denotes the elementary charge. Expanding the various terms in the effective Lagrangian in powers of the pion 
field and using the canonical formalism along the lines of Ref. (47j , we end up with the following interaction terms in 
the Hamilton density 

n<§ = ^ (rf [f x fr] N) ■ (iVt [f x fr] N) , (2.15) 
and the electromagnetic current density is of the form 

J°J~ 1] = |iVt(l + T3)JV, 

J^- 1] = e[ 7 ?x#] 3 , 
J02 ( 1} = -e[n x Vtt] 3 , 
J 2 i (0) = e^N^a[TX7f] 3 N. (2.16) 

In the above expressions we adopt the notation of Ref. [3^|. In particular, the subscripts a and b in TL^ and J^ b 
refer to the number of the nucleon and pion fields, respectively, while the superscript k gives the dimension of the 
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Class 1: 



Class 2: 



Class 3: 



Class 4: 



Class 5: 



Class 6: 



Class 7: 




FIG. 1: Diagrams showing contributions to the leading two-pion exchange currents. Solid and dashed lines refer to nucleons and 
pions, respectively. Solid dots are the lowest-order vertices from the effective Lagrangian while the crosses represent insertions 
of the electromagnetic vertices as explained in the text. Diagrams resulting from interchanging the nucleon lines are not shown. 



operator as defined in Eq. (|2.6[) . Further, the symbol 
spaces. 



in Eq. (|2.15p denotes a scalar product in the spin and isospin 



The formal operator structure of the leading two-pion exchange two-nucleon current at order 0(eQ) is given in 
appendix [5] For the sake of convenience, we distinguish between seven classes of contributions according to the 
power of the LEC <?a (i.e. proportional to g\, g\ and g\) and the type of the hadronic current J^ , or Jq 2 as 
shown in Fig. [T] Notice that there are no contributions proportional to g° A and involving J 20 and J 21 . We also 
emphasize that the second diagram in the class 3 does not generate any contribution. It results from the term in the 
Hamilton density which is absent in the Lagrangian and arises through the application of the canonical formalism. 
Finally, it should be understood that the meaning of diagrams in the method of unitary transformation is different 
from the one arising in the context of covariant and/or time-ordered perturbation theory. The diagrams shown in 
Fig. [1] serve merely to visualize the topology corresponding to a given sequence of operators H and J appearing in 
the formal expressions given in appendix [SJ 

It is a fairly straightforward albeit tedious exercise to evaluate the contributions to the nuclear current corresponding 
to the operators given in appendixfSJ Below, we give explicit results for the current and charge densities, — (p, J), 
resulting from the individual classes using the notation 



(pi'pbVHPi^) = 5{Pi +P2 -Pi -P2-k) 



' c7 

£ 

.X=cl 



J£ + (1 «-> 2) 



(2.17) 
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where ft (ft') denotes the incoming (outgoing) momentum of nucleon i and fc is the photon momentum. Further, 
(1 <-> 2) refers to the contribution resulting from the interchange of the nucleon labels. We find the following results 
for the current density from the individual classes of the diagrams shown in Fig. [TJ 



Jcl (51,32) 



ft2(ft,ft) = 



9aJ_ 
16Fi 



ft [n x f 2 ] 3 + 2 [ft x ft] t( 



cri 



9A i 

647r 2 Fj 



ft [n x f 2 ] 3 + 2 [ft x ft] rf 



(27r) 3 + w_) 



4 ■ 

9a* 



<i 3 Z cftl + cj + o>_ + cj 2 



4r| / (I • [ft x ft] ) - (Z 2 - q 2 ) (ft [ft x f 2 ] 3 + 2 [ft x ft] rj 3 



16F4 y (2tt) 3 w 3 u 3 (w+ + w_) 

ft [ft x ft] 3 (8M 2 + 3g 2 ) + 2 [ft x ft] t 3 (8M 2 + 3g 2 ; 



4 ■ 



64tt 2 F4( 9 2 + 4M2) 



4 (qf + 4M 2 ) [<h x *i] 



Jc3 (ft, 92) 

Xl (ft, ft) 
ft 5 (91,92 ) 

^c7 (91,92) 



in x r 2 



128F4 

x (k 2 + fa) , 

Jc6 (ft, ft) = 0, 
g\i f d 3 l 



d 3 l 



L(qi), 

1 / (uJi — W 2 ) (Ml — (J3) (u)\ + UJ 2 )((jJl — W3) (Wl — W2)(W1 +W3) 



(27l) 3 WiW 2 W 3 \ (wi + U)2)(u)l + Cl> 3 ) (wi + w 3 )(w 2 + W 3 ) (W1+W 2 )(W2+W 3 ) 



CJl + UJ 2 + W3 



y (27r) 3 wiw 2 w 3 (cjj + w 2 )(ui + w 3 )(w 2 + w 3 ) 



(ha + k 3 ) 



[ft x t 2 ] 3 ft • fc 3 



2r 3 ft 



4 • 



ft X fc 3 

d 3 z 



lJ 2 



y (27r) 3 wi^ 2 w 3 \wio; 2 cj 2 w 2 ft> 2 — + lo 3 ) cj| (oft - wf) (^2 + ^3 ) 



x ( fc 2 + 



[ft x ft] fci • fc 2 fci • fc 3 + 4rf ft. 



fci x fc 2 



fci • fc 3 



(2.18) 



while the contributions to the exchange charge density read 
Pel (ft, ft) 
Pd (91,92) 



Pc2 (ft, ft) = Pc3 (ft, ft) = , 
d 3 l 1 



e^-r 3 



(2tt) 3 w 2 w 2 



[J 2 ~ 2 



9iJ 



327rFi 



r 3 [M T + (2M 2 + g 2 ) A( qi )] , 



Pc5(ft,ft) = 



Pc6 (91,92) 



-eAr 3 



16F4 

32^ 

9 a 



d 3 l 



(2tt) 3 w 2 cj 2 



V 2 



r 3 [M w + (2M 2 + g 2 )A( g i)] , 



d 3 J 



8F*J (2tt) 3 



At 3 1 • [ft x ft] T [ft x ft] + (l 2 - qf) t 3 



9% 



327r.Fi 



9 A 



#=7(91,92) = e— j 



(4M 2 + 2g 2 ) t 3 + r| [ft • ft ft • ft - ft* ft • ft] ) 
d 3 Z 1 



8F4 y (2tt) 3 oJiJiJi 

[ft x f 2 ] 3 fci • fc 2 fci 



( T i + T i) [ki ■ k 2ki ■ fc 3 + fci 

fc 3 X ft 



fc 2 x ft 



fci 



(llM 2 + 3g 2 ; 
4Af2 + g 2 

fc 3 x ft 



(2.19) 



where ft = pi' — pi is the momentum transfer of nucleon i. Clearly, ft and ft refer to the spin and isospin Pauli 
matrices of the nucleon i. Here and in what follows, we label the Cartesian components of various vectors in isospin 
space by the superscripts rather than subscripts in order to avoid a possible confusion with the nucleon labels. Further, 
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lu± = \J (I + <fi) 2 + 4M | and = \Jkf + M 2 with fcj = fc 2 = l — q\ and fc 3 = Z + 9V We also use here the short-hand 
notation with qt = |^|. Finally, it should be understood that the above expressions refer to matrix elements with 
respect to momenta and operators in spin and isospin spaces, see Eq. (|2.17|) . The expressions for J c i, C 2 (91,92) and 
Pc4,c5,c6 (91, 92) only involve two-point functions so that the corresponding integrals can be evaluated in terms of the 
loop functions L(q) and A(q) defined as 

Is / s + q 



L(q) = --In ^ , withs= vV + 4M 2 , 
z q \s q J 

A(q) = ^arctan(^-). (2.20) 

Notice that the loop integrals are ultraviolet divergent and thus need to be renormalized by an appropriate redefinition 
of the LECs accompanying the NN contact interactions. Here we only show the non-polynomial parts of the resulting 
expressions which give rise to long-range contributions and are uniquely defined once the regulator (i.e. cutoff) is 
removed. Short-range contributions involving NN contact interactions will be considered elsewhere. The symbol — > in 
the above equations signifies that the original ill-defined expression is evaluated using dimensional regularization 1 and 
only non-polynomial in momenta and M 2 contributions are shown explicitly. The exchange current/charge density 
contributions J C 3, C 5,c7 (9*1, 9*2) and p c r (q\, 92) result from loop diagrams where a photon couples to pions in flight. The 
corresponding loop integrals depend explicitly on two external momenta 9*1 and 92 and can be written in terms of the 
three-point functions. The complete expressions in momentum space are given in appendix [B] in terms of a number of 
scalar integrals, which are easily calculable numerically as detailed in appendix [C] We also emphasize that we have 
explicitly verified in appendix |D] that the derived exchange currents fulfill the continuity equation. Notice further that 
all loop integrals can be carried out analytically in configuration space as discussed in the next section. The results 
for the current density presented here agree with the ones of [34|. 2 



III. TWO-PION EXCHANGE CURRENT IN CONFIGURATION SPACE 



The obtained expressions in momentum space depend only on the momentum transfers of the individual nuclcons 
leading to a local form of the current operator in configuration space which is defined according to 

J»(no,r 20 ) ~ J e*-*>J"(fi,fi), (3.1) 

where r 10 = Pi — r , r 20 = ^1 ~ ^0 an d f\ , r 2 and r denote the positions of the nucleons 1 and 2 and the photon coupling, 
respectively. Using the formulae collected in appendix [E] and expressions for the current operator in momentum space 
in terms of three-dimensional loop integrals given in the previous section, we obtain the following surprisingly compact 
expressions for the current density 



X2 (fio,f 2 o) 

J C 3 (no, ^20) 



128tt 3 F4 
„ 9 a Ml 



V10 [fi x f 2 f + 2 



V10 X <7 2 



iu(2xi ) 



32^ r 10 x 

Ml 
5f2vr 4 J F 1 i 



V10 [fi x f 2 ] J + 2 

0(X20) 2 I 



V10 x a 2 



S(x 2 o) 



K Q (2 Xl 



r ^ -.13 /rf y4 \ K 2 (x W + X20 + X12) 

[n x r 2 ] (V10 - v 20 )- 



(Xio X 2 Xi2)(x W + X20 + Xu) ' 



1 Using e.g. a cutoff regularization will lead to the same result after taking the limit A — > 00. 

2 For contributions involving the three-point function we could only check the intermediate results for the integrands involving different 
pion energies. 
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J C 5 (no, r 20 ) 



Jd (rio,r 2 o) 



9\Ml 



Xi(xio + x 2 o + X12) 



[ri x f 2 ] V12 • V20 - 2Tf<7 2 



V12 x V 



20 



(Xio X20X12) 

^(V 10 -V 2 o 



512 



[n x f 2 ] 3 V12 • Vio V12 • V 20 + 4t 2 3 cti 



V12 x V 



10 



V12 • V 



2(1 



£10 + 2>20 + Zl2 

x Ao(a;io + 2:20 + £12) 

2^10 ^20 ^12 



(3.2) 



and the charge density 



Pel (n , f 2 o) = p C 2 (no, r 2 a) = Pc3 (rw,r 20 ) = , 
Pc4 (r 10,^20) 



e ji^rf^o) (V? -2) 



2567r 2 Fi 



'10 



Pes (rio,^2o) 

Pee (no, ^20) 
Pc7 (no,r 20 ) 



t 2 d(a; 2 o) ( V 10 



256tt 2 F4 
256tt 2 F4 
128tt 2 F4 
5127r 3 Fi 
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5 (£20) 



if (2V 2 



^20)^ (3V 2 -11) 



r| (?i ■ Vio (?2 ■ Vio 

,-2x 10 



r 2 a± ■ <J2 



-2x x 



10 



no 



(Tf+T 3 ) Vi2-Vi Vi2-V 



20 



[fi x f 2 ] V12 • Vio V12 



V20 x 02 



V12 • 











Vio x cti 


V12 • 


V20 x a 2 


) 



Xw X20 X12 



(3.3) 



In the above expressions, K n (x) denote the modified Bessel functions of the second kind and we have introduced 
dimensionless variables x\q — M^fio, X20 — M^^o and X12 = M^f^ — M^ifx — ). Further, x%j = \xij\ and all 
derivatives with respect to xio, X20 and X12 are to be evaluated as if these variables were independent of each other. 
We also emphasize that the above expressions are valid for xio + £20 > 0. Finally, it should be understood that the 
behavior of the current and charge densities at short distances will be affected if one uses a regularization with a finite 
value of the cutoff. 



IV. SUMMARY AND OUTLOOK 



In this paper, we applied the method of unitary transformation to derive the leading two-pion exchange two-nucleon 
charge and current densities based on chiral effective field theory. The resulting nuclear current is given both in 
momentum and configuration space. The results in momentum space involve the standard loop functions L(q) and 
A{q) and, for certain classes of diagrams, also the three-point functions. In the latter case we exp ressed all tensor 
integrals in terms of a set of scalar integrals in different dimensions similar to the method of Ref. [4a |. see appendix [Bl 
For the remaining scalar integrals we derive in appendix [C] a simple representation using Feynman parameters which is 
well suited for numerical calculations. In configuration space, we were able to evaluate all loop integrals analytically 
leading to very compact expressions in terms of the modified Bessel functions of the second kind. We have also 
explicitly verified that the derived exchange currents fulfill the continuity equation, see appendix ID] 

In addition to the two-pion exchange contributions, there are also one-pion exchange and short-range terms at order 
O(eQ), see Refs. (33l . l34j for a recent work based on time-ordered perturbation theory. As will be demonstrated in 
a subsequent publication [46| , renormalization of the one-pion exchange contributions at the one- loop level strongly 
restricts the ambiguity in the definition of the current and provides a highly non-trivial consistency check of the 
calculation. In particular, one needs to ensure that all appearing ultraviolet divergences are absorbed into redefinition 
of the LECs di and k from C^ N and respectively, with already known /3-functions, see e.g. [4!| [5(| HH, ■ This 
work is in progress [4q |. 

Finally, in the future, one also needs to test the convergence of the chiral expansion for the one- and two-pion exchange 
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currents by calculating the corrections at order 0(eQ 2 ). Given the large numerical values of the LECs 03,4 from C^h, 
one might expect sizeable corrections which, indeed, is well known to be the case for the two-pion exchange potential 
[53l ]. In this context, it might be advantageous to include the A(1232) isobar as an explicit degree of freedom in 
effective field theory utilizing the small scale expansion [54[ , see [5^, [Hfl, [57], [H, [5!^ [6(| for recent work along these 
lines in the purely strong few-nucleon sector. 
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APPENDIX A: FORMAL STRUCTURE OF THE LEADING TWO-PION EXCHANGE CURRENT 

OPERATOR 



In this appendix we give the formal structure of the two-pion exchange current operator that results after applying 
the unitary transformation. All vertices entering the expressions below should be understood as second-quantized 
normally-ordered operators. Here and in what follows, we adopt the notation of Ref. [39l | with a few minor modifica- 
tions. 
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APPENDIX B: LEADING TWO-PION EXCHANGE CURRENT IN MOMENTUM SPACE 



In this appendix we give the expressions for the two-pion exchange current operator in momentum space. Following 
Ref. [6lj . the most general expression for the current and the charge density can be written as 



5 24 
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where f- = f- (gl, 92) are scalar functions and the spin-momentum operators Oi and Of are given by 
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Of = 1, 

Of = 5*1 • [52 x o 2 ] + 5*1 • [5*2 x oi] , 

Of = 5*1 • [52 x ct 2 ] - gl • [g 2 x cti] , 

Of = 01 • 02 , 

Of = (51 • o* 2 ) (52 • 01) , 

Of = (51 • 0*1 ) (52 • o 2 ) , 

Of = (52 • 01) (g 2 • o 2 ) + (51 • 01) (gi • o 2 ) , 

Of = (52 • 01) (g 2 • o 2 ) - (51 • 01) (gi • o 2 ) • (B.3) 



As a basis for the isospin operators we choose 
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The nonvanishing long-range contributions to the scalar functions f- = f- (q\, q 2 ) due to two-pion exchange calculated 
using dimensional rcgularization are given by 
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In addition, there are nonvanishing functions f 2 given by 
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In the above equations, z = q\ ■ q 2 , qi = \qi\ and the loop functions L(q) and A(q) are defined in Eq. (|2.20|) . Further, 
the functions / correspond to the three-point functions via 
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Here, all propagators are understood to have an infinitesimal positive imaginary part. Notice that here and in what 
follows, we will only need the functions I for four-momenta with vanishing zeroth component. We further emphasize 
that all functions l\^^ 2 „ 3 ) which enter the above equations except I^f^X, /(2 To) anc ^ {120) arc nn it c m dimensional 
regularization in the limit d — *• 4. For these functions, we define reduced functions by subtracting the poles in four 
dimensions 
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Here, /1 is the scale introduced in dimensional regularization and 7e = — r'(l) ~ 0.577. Finally, for scalar functions 
contributing to the charge density we obtain the following expressions: 
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APPENDIX C: EVALUATION OF THE 3-POINT FUNCTION 



The momentum-space expressions given in the previous appendix involve the 3-point function defined in Eq. (|B.8j) . 
Below we show how loop integrals of this kind can be evaluated by introducing the corresponding Feynman parameters. 
In particular, we consider the following integrals 

Ei = 7(d;0,^i;p2, ^25^3,^3; 0,^4), 

I 2 = l(d;0,vi;p 2) v 2 \V^z\ 0,0), (C.l) 

in the following kinematics: 

qi = (0, qi), qx ■ qi = z, v 2 = 1, V ■ qi = v ■ q 2 = 0. (C.2) 

Here, we consider the case with ^ = 1,2,3,... and | > 0. Starting from this point, we denote by % the length of the 
corresponding three-momentum in units of the pion mass, i.e. qi = \<ji \/M K . Introducing the Feynman parameters 
and carrying out the integration over the loop momentum, we obtain the following result for the first integral in 
Eq. JCU): 
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where 



and we have introduced 



E E 2 + A AD E E 2 + A AD _ . 

yi= 2D + y AD 2 ' V2= 2D-i AD 2 ' (C ' 4) 



A = l + q((l-t)t > 0, 
-S = 2qiq 2 z + q2, 
C = -2q x q 2 z, 
D = q\ > 0, 

£ = B + tC. (C.5) 
It can be shown that the following inequalities hold true in the integration region for the considered kinematics: 

yi >t > y > > y 2 . (C.6) 
Thus, the remaining two-dimensional integral in Eq. (|C.3P can be easily calculated numerically for all desired values 

Of Vi. 



Similarly, for the second integral in Eq. (|C.1|) we obtain: 
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where, again, the inequalities given in Eq. (|C.6[) hold true. 

For the reduced functions ^^(1 1 o)' ^rccT (1 2 0) anc ^ ^rccT (2 1 0) we obtain the expressions: 
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APPENDIX D: CURRENT CONSERVATION AND THE CONTINUITY EQUATION 

Current conservation implies that the electromagnetic current operator J(x) should fulfill the continuity equation 

V • J(r) = = -i [H, p] = -1 [H +V,p], (D.l) 

where p — p(f f ) = J°(r) is the charge density and H, Hq and V refer to the two-nucleon Hamilton operator, kinetic 
energy term and the potential, respectively. The continuity equation thus provides a powerful check for the calculation. 
For the leading two-pion exchange contributions to the current operator the continuity equation takes the form 

V ■ j£\r) = e [n x f 2 f W^{?x - r 2 ) [S(r- n) - 5(f- r 2 )} , (D.2) 

where the potential in the isospin limit is written in the form 

V = V + fx ■ f 2 W , (D.3) 

r, denotes the position of the nucleon i and the superscripts of J and W refer to the chiral order. Further, we made 
use of the explicit form of the leading one-body charge density 

P(r) = \ [(1 + rf) 6(f- n) + (1 + r 2 3 ) S(f- r 2 )] , (D.4) 

and the absence of the leading one-pion exchange charge density, p\^ = 0. Eq. (|D.2[) can be transformed into 
momentum space 



i(qi + &) ■ J2\<h>&) = e [n X f 2 ] 6 \W%>(&) - W£\q 2 )\ . (D.5) 

The left-hand side of the above equation can be expressed in terms of the basis operators Of defined in Appendix [B] 
Using the representation of the current operator in Eq. (|B.1[) , the spin- momentum operators appearing on the left-hand 
side of the above equation can be expressed in terms of the operators Of as follows: 

24 8 

SI (quQ2) k-6 3 = 5>f («7l,g 2 ) Of , (D.6) 
where the scalar functions gf S (<fi, q 2 ) read 

9f = -ft ~f!+ \k 2 fj + \(ql- ql)ff + \k 2 ff -\{q\- ql)f}\ 



(2)^ w( 2 ) 
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gf = -f! +f! + \k 2 f! + - q M + \k 2 fl° - - <i)fl 
gf = *V, U + 

rf* = ^ 17 + (^-gI)// 8 + 2/f 1 , 

ft 75 = fl 9 + /f + ~* 2 // 3 + 5(9? - 9 2 2 )// 4 + ^ 2 /P + - ll)f?\ 

gf = -f 20 + /f - ~fc 2 // 4 - ^fe 2 - ? 2 2 )/ 13 + ^ 2 / 24 + \{dl - al)f 23 , (d.7) 

and k = <fi +q 2 . Substituting the expressions for ff for the leading two-pion exchange current operator from Eq. (|B.5[) . 

iS (2tt) 

one finds that the only non-vanishing function g { is given by 

^ = e 3^{ 4M ' W ~ ^ -^ +q * ^ ~ 10 & - x ) + wtimz} +a ~ (1 ~ 2) • (D - 8) 

where a denotes a polynomial in qi and 92 whose form depends on the choice of the subtraction scheme. Using the 

(2) 

explicit form of IU 2 V , 



w?J = - L{q) 



2w 3847r 2 Fi 



{iMl (5g\ - 4g 2 A l) + , 2 (23^ - 1 £ - l) + §§f^) + , (D.9) 



where a' again denotes a subtraction-scheme dependent polynomial contribution, it is easy to see that Eq. (|D.5p is 
indeed fulfilled. 



APPENDIX E: COORDINATE-SPACE REPRESENTATION OF VARIOUS LOOP INTEGRALS 



In this appendix we collect the formulae needed to obtain the expressions for various loop integrals in coordinate-space. 
The two-pion exchange contributions to the current operator resulting in the method of unitary transformation are 
given in terms of three-dimensional loop integrals with the integrands being rational functions of the pion energies 
u>i = \J kf + A/ 2 where ki = \ki \ refers to the pion momentum. In order to end up with simple expressions in 
coordinate space, it is convenient to rewrite the integrands as continuous superpositions of the propagators using the 
following expressions: 

1 2 f°° 1 1 

d(3- 



^1^2(^1+^2) ir Jo ujf + (i 2 uol + P 2 ' 

U)\ + UJ1LO2 + U>2 Id 1 

lo\u^(uji + 012) 2 <9M 2 ^1^2(^1 + 0J2) ' 

= *l d(3 {(^ + m4 + p 2 ) + R + p 2 )^i + p 2 ) 2 1 ' (K1) 

In addition, we also need the following relations which involve three different pion energies: 

1 UJl + UJ 2 + ^3 



D x = = 



UO1UJ2OJ3 (UJI + UJ 2 )(UJ1 + W 3 )(lJ2 + W3) 



2 f°° 1 

d/3 



Wo -( w 2 +/32)(^+/32)(^ + /32)' 
D _ 1 / (Ul - ^2X^1 - ^3) (^1 + ^2X^1 - ^3) (^1 -^2X^1 +^3) 

W1W2W3 \ (Wl + W2)(wi + W3) (U0\ + UJz)(uJ2 + W3) (ll>1 + W 2 )(w 2 + UJ3) 

4 [°° ( 4/3 2 1 

dp 



TV 



D* = 



(coi + p 2 m + o 2 m + p 2 ) k 2 + Pm + p 2 ) 

1/1 U>2 OJl 



LO1LO2UJ3 w 2 (oj 2 — + cj 3 ) w|(u; 2 — w|)(w2 + W3) 
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dM%0JiLl2W3 [U3\ + u 2 )(oj 1 + L0 3 )(u0 2 +W3) 



rf/3 



1 1 

+ 



(a, 2 +/3 2 ) 2 ( W 2 + /? 2 )(^ 2 +/3 2 ) (u; 2 + /}2) ( w | + £2)2 (^2 + £2) ( W 2 + 02)( w 2 + £2)^2 + £2)2 

Using the above expressions, it is straightforward to carry out the Fourier transformation !F, 

T {m , q 2 )),/ ^|^e*-e- - /(* , ft ) , (E.3) 
of the integrals which appear in the calculation. For example, 



1L I ^ = s(f 20 ) f —tL [*L e **»l r d 8 I 

(2^ w+w _( W++W _)j 1 20) J (2tt) 3 J (2tt) 3 irJo P (cu 2 + + /3 2 )(a, 2 + /3 2 ) 







M 3 roc e -^/i+^2xi 

= / d ^ — sis 



(2tt) 3 (2tt) 3 (4w 2 + /3 2 )(4w 2 + 

e~ 



4tt 3 y 2x 2 

- s^^ 20 )^^' (K4) 

where cj± = \J (I ± q\) 2 + 4M 2 and the i-Q refer to the modified Bessel functions. We made use of Eq. (|E.1[) in the 

first line and carried out the substitutions t + q\ — > 2fci, £ — <fi — > 2fc 2 in the second line. Notice further that the last 
equality is valid for x\o > 0. In a similar way one obtains 

f f dH I 2 \ Ml 2 K&x w ) 

V7 (2?r) 3 w 3 w 3 _(cj + / 16?r 3 v ; xio 

Again, these expressions are only valid for xw > 0. Clearly, the above results may also be obtained by first carrying 
out the momentum-space loop integrals. The resulting non-polynomial expressions are given in terms of the loop 
functions L and A, see Eq. (|2.20[) , and can be easily Fourier-transformed using the following relations 



(2tt) 3 s 2 4tt x ' 



j3 M" 2 +2A; p- 2:E 

' 9 e*V)M(g) - ( -l)fc(V 2 -4) fc i^V' (E-6) 



(2tt) 3 v ' v / v x , 2ii 

(2tt) 3 ~ v ~ ' "™ K y ' x 8vr x 2 

where k = 0, 1, 2, 

Consider now the more complicated expressions involving three different pion energies which appear in the calculation 
T \ I T^^^fe S(k 2 -h+ q\) 6(k 3 -hi- q 2 ) (k 2 + k 3 ) D 2 
d f3 I f\f\ d \ e -^-r w ^ 0e 



(2tt) 3 

4 r jo f d3k i d3k 2 d 3 k 3 __ig 2 . rl „ ^ ^. ? J 4/3 2 J 



7T 

x (fc 2 + k 3 ) 



(2vt) 3 (2,t) 3 (2tt) 3 V(^ 2 +/3 2 )(^ 2 2 + /3 2 )K 2 +/3 2 ) (w 2 2 +/3 2 )K 2 +/3 2 ) 

j_/Vf7 roo g-Vi+Z^zio e -V 1 +/ 32a; 20 e -Vl+/3 2 ^i2 



»7 /-OO 

f (Vio- V 20 ) y d/3/3 : 



4tt 4 Jo scio x 20 X12 

iMj - ^(xiq +x 20 +x 12 ) 

-— (Vio-V 20 )7 T7 ■ ■ r, (E.7J 

4tt 4 {x w x 2 oXi 2 ){x w + x 20 + xi 2 ) 
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where in the first line we used Eq. (|E.2|) . The last relation is valid for positive values of the argument of the Bessel 
function. Similarly, we obtain: 



fe, —Afc-KlfalO + ^20 +X12) 



M 4 (iM \ n + m +P 

V J Q0^4 V 12 V 12 v 10 V 10 V 20 V 20 /_ _ _ \ ' 

(a;i0^20^12) 

J 7 (/ J^ d3k2d3k3 S ^ 2 ~ kl+ ^) _ * l _ « 2 ) ^ ' ' ' fc l Wfc 2 1 ' ' ' ' ' ' k 3 P D 3 

(-1) +p V 1 1 2 ---V 1 ^Vf ---V J 1 'oV 2 ^---V 2 g ^0(210 +X20 +X12). (E. 

6ZTT^ XlQ X20 X\2 
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